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TECHNICAL MEMORANDUM 

A PRELIMINARY LOOK AT CONTROL AUGMENTED DYNAMIC 
RESPONSE OF STRUCTURES 

I. INTRODUCTION 


Classically, structural dynamieists and structural engineers have looked at structural redesign as 
the technique for solving structural response and loads problems. Generally, control engineers have used 
their discipline as a means of controlling some system state, such as space vehicle attitudes, automobile 
pollution control, etc, In recent years, the two fields have started merging when control has been used 
as a means of reducing structural weight in the aeronautics, aerospace, and automobile industries, The 
Space Shuttle is a good example of the uso of vehicle load relief control system approaches to reduce 
the overall vehicle aerodynamic loading and thus save structural weight, Tins system used pitch and yaw 
acceleration feedback to reduce aerodynamic loading by reducing angle of attack, side slip, and rolling 
of the vehicle in such a manner as to load the Orbiter wings in the most favorable direction. In addi- 
tion, elevon load relief was employed to reduce elevon loading during ascent, This system, however, 
only treated vehicle rigid body response. Improved approaches are being developed which move beyond 
rigid body response dealing with elastic body response leading to control configured, optimized design 
configurations. Active control configured flutter suppression, aeroelastic tailoring, modal suppression, 
and optimal design techniques fall under this general heading. These techniques have evolved m struc- 
tural dynamieists and control engineers have recognized the potential of using control systems and 
control logic as means of altering structural dynamic responses, thus replacing structural redesign or 
structural weight with more efficient use of already existing control systems. Optimized design of new 
configurations naturally follows through the use of more complex control systems. 

The evolution of the multidiscipline, structural control interaction has not developed as fast as 
it should. Several reasons for this slow development are clear: (1) normal protection of one’s own 
discipline, (2) the uso of different transformations, terms, nomenclature, etc., for solving the same 
type differential equations, and (3) lack of proper systems engineering and organizations to force the 
cross fertilization and system trades across these major disciplines. 

Classically, if a structural dynamicist wants to change response or reduce loads, he changes stiff- 
ness by (1) adding or subtracting materials, (2) passive isolation of components, etc., (3) addition of 
passive dampers, and (4) detuning the system from the forcing function, This report will treat these 
same concepts first from classical vibration and elasticity theory and then show how control logic can 
accomplish the same goals while still preserving the nomenclature and form structural dynamieists are 
familiar with. A brief look will also be made in the control engineer's field showing how, with mini- 
mum effort, one can transform the knowledge and insight from structures to control and vice versa. 

This is accomplished by first repeating the basic vibration characteristics of a single degree-of-freedom 
mass, spring, damped, forced system, The basic equation is then recast with a control system feedback 
logic, then put in the same form as the original equation, thus preserving the response characteristics 
a structural engineer is used to working with. In the reformulation with control logic, the basic parame- 
ters of inertia (mass), damping, and stiffness are augmented with control parameters. Next a single 
bending mode is treated to show the transference of the single degree-of-freedom to include structural 
gains, central control, and distributive control concepts. Finally, these same concepts are briefly looked 
at for two- and multimode-structurul systems. 


ii. SINGLE DEGREE-OF-FREEDOM VIBRATION SYSTEM 


The linear single dcgrec-of-frecdom forced vibration system has been analyzed, discussed, and 
published probably more than any other dynamic system. The reason for this is obvious; it is the basis 
for dynamics and vibration in general. The same is true for linear proportional gain control theory, 
Since the purpose of this report is to put fundamental control theory in the vernacular of the struc- 
tural dynamicist, a brief review of the basic results of vibration theory will be stated. Additional 
information can be found in any vibration, dynamos, or control theory text book, Figure l depicts 
the classical single degree-of-freedom linear forced system. Constant properties are assumed. 



Figure 1. Single degrco-of-freedom system,, 
The describing differential equation is; 


mx + cx + Kx - f(t) 


( 1 ) 


Solutions to linear differential equations in this form break out into two parts, homogeneous and 
nonhomogenous. Out of these have grown the classical stability criteria, usually expressed in structural 
dynamics as frequency and damping and in control theory as complex roots, the real part depicting a 
product of damping and frequency and the imaginary part as the damped frequency. Without repeating 
all the derivations, equation (l) can be rewritten in the form; 


x + 2 frox + cA = f(t) (2) 

n n 

where, 

\ f(t) C 0 IC 

t(t) = — and 2 fco [ - — <of - — , (3) 

m 11 in 11 m 



and in the homogeneous form where f(t) - 0, the classical roots are 6 + ice, where, 


5®-£co 


a>d ~ cu 




(4) 


Solutions to this equation for various type forcing functions are well known and are familiar 
to the control and structural dynamics community. 

The solution for forced oscillation can be obtained in several forms. The results are presented 
as classical solution curves parameterized in terms of damping and frequency. In general, these solu- 
tions are expressed in terms of forcing functions that have the form of (1) impulses, (2) steps, (3) 
sinusoidal, (4) ramps, etc. The response, for example, to a step function, subcritically damped, is well 
known and takes the form: 


X(t) « x st 




(5) 


where X s t is the new equilibrium position and 

£ 

tan a = (6) 


The maximum displacement is then 


^\max “ ^st 




1 + — cos a 

%/ 1 -? 2 


(7) 


Plotting the solution for various damping values is shown on Figure 2. 

It can be seen for increasing damping, the subcritical damping £ < 1 , the extreme values shift to 
the right (frequency decreases), and the amplitude decreases. If f(t) is a harmonic forcing function, the 
homogenous solution takes the form 


X(t) 


(V 0 + & x 0 y 


x° + 


oj 2 (1 - £ 2 ) 


e -£a>t cos / w y/ 1 _ 


It + a 


( 8 ) 


Figure 2, 


Displacement of a subcritieally damped vibration system to step input, 

,vith 

-(V 0 + X 0 M (9 

tan $ a — . — ~~*2=ssssss=i 

Xococ/r? 

The solution in graphical form is shown in Figure 3, 



The nonhomogenous solution is better approached from the vectorial form if harmonic forcing 
functions are assumed. 

MX + CX + KX~F 0 e int 


( 


Thus, the particular solution is of the form 


X*X e int 


( 11 ) 


giving 


x [-a 2 + me + k) « F 0 


( 12 ) 


EquationJlB) shows the equilibrium of the various forces, inertial, damping, restoring, and driving, 
with the F 0 and X* being no longer collinear, as in undamped oscillations, but forming a phase angle, a. 

Letting r 3 £2/co be the ratio of the forcing frequency to the undamplcd natural frequency co, the 
equation becomes 




iy K 

“ — 

<*■> J, O T 

- r") + 4 f' " r“ 



(13) 


and 


C S2 

tan o: ~ - 

k - m a- 


Plotting ju, the classical magnification factor gives the valuo as shown in Figure 4 and the phase angle 
as shown in Figure 5. The peak magnification factor ju max and /u versus damping ratio is shown in 

Figure 6. From these curves, one can sec the classical changes in response familiar to ail structural 
dynomicists and control engineers, 



Figure 4. Magnification factor. 



c 



Figure 5, Phase angle ety 
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Figure 6, Damping factor 


The effects of dumping, inertia, and stiffness upon the response as well as how they interact 
with the forcing function are well known, For example, more damping reduces the response but 
increases the period of the oscillations, more mass reduces the static responses but increases the period, 
increasing the spring increases the frequency, and decreasing the period shortens the deeay time, 

Finally, tuning the frequency of dynamic systems to the forcing function increases several fold the 
response amplitude, The structural engineer by changing mass, damping, and stiffness alters the response 
of the system, Isolation systems are derived using these concepts, In summary, the response can be 
changed by (l) adding mass, (2) adding structural damping, (3) changing stiffness through material 
changes, structural configuration changes, etc,, or (4) altering the forcing function or detuning the struc- 
tural response from it, 


Section III shows how the system response alteration can be achieved by augmenting it with 
control concepts instead of structural changes, 



b 


SECTION III. SINGLE DEGREE-OF-FR c .EPOM VIBRATION SYSTEM WITH CONTROL 


From control theory, forces can be generated which are linearly proportional to the displacement, 
velocity, and acceleration. Using the familiar SDOF, tho application of control theory to alter system 
response is illustrated in Figure 7 where 


5(t) - l(a 0 x(t) + a 1 x(t) + a 2 x(t)l F* 


(14) 




¥ 




Figure 7, Mass spring system with control feedback. 


A control force is generated linearly proportional to displacement, velocity, and acceleration, and 
f(t) is a completely independent forcing function (non-feedback) of the system. For tliis formulation, the 
assumption is made that the control force is ideal; i.e,, no phase lags exist in the control mechanism. The 
introduction of phase lag or lead can be added later to further illustrate control augmentation of dynamic 
response of structures, Introducing the augmented control force into the equation of motion of 
dynamical systems gives: 


nix + cx + Kx - f(t) + 5 (t) - mx + cx + Kx « f(t) + (oqX + apx + a 2 x) F* 


(15) 


Recombining or collecting terms, 


(m - a 2 F’ ,< ) x + (c - a^F 11 *) x + (K - agF*) x = f(t) 


(16) 


It is clear that equation (16) is identical in form to equation (2). What has changed is the definition of 
the parameters. No longer Is the mass term true mass, but it is an apparent mass due to the control 
force augmentation. The same is true for damping and stiffness terms. In addition, tho control parame- 
ters (gains) can take on positive or negative values, hence they must be carefully chosen in order to keep 
the system stable since it is possible to drive the augmented terms to negative values. Defining the 
augmented terms as: 
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( 17 ) 


m* « ni - ajF* 
c* » c ~ ajF* 
K* «= K - a 0 F* 


equation (16) becomes: 


m* x + c* x + K* x 13 f(t) 


08) 


where, 


o a -2fw 



c-ajF* 
m - aJF* 


|f H 1 K - UqF* 
in* m - a 2 F* 


09) 


The ability of dynamic engineers, through application of control augmentation, to alter the 
behavior of dynamical systems while presewing the structural design and geometry is clearly illustrated, 

In general, control theory and control systems are thought of as tools for achieving a desired 
response of a system already in design, verification, or operational phases, These uses of control employ 
both active feedback and open loop command control techniques. The past few years, systems engineers 
have recognized the power of control to achieve the response goals and supplant in a more optimum way 
structural weight, stiffness, and damping. Tills allows for tighter design tolerances and higher perform- 
ance systems at lower cost and risks, This is accomplished through the use of the extra variables intro- 
duced by control, including the choice of locating forces optimally for desired response, This will be 
clearer in the sections dealing with space vehicle elastic modes and non-ideal control. This approach has 
been exemplified extensively in the active flutter suppression techniques employed on modern aircraft 
and the new field commonly called “control configured design,” What is being done essentially is using 
control theory prudently to augment structural design parameters, such as stiffness, damping, and mass 
in both a static and dynamic sense and detuning the system from the forcing function. 

The example given used only one sensor of each type and one control force; however, for a multi- 
mass system, one can go to multisensor, multicontrol force systems and extend substantially the number 
of variables for optimizing the system. In the past, structural response has been augmented using a 
control system designed for another function; orientation of space vehicle along a desired flight path. 
Using the concepts presented, two control systems or an augmented single control system can be made 
more optimum to perform the two separate tasks, orienting a vehicle’s attitude and changing the struc- 
tural dynamic response to achieve a more optimum structural design. 
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The substitution of equation (10) is not made back into the equations, since it is obvious how it 
is accomplished without destroying the validity of the solutions or graphical solutions shown in Section 
II, 


At this point, it is clear that a control system can be used to act ns structural elements and alter 
the system responses, Classically, many ways are available for dealing with and understanding the system 
response and determining the best approach for arriving at a solution, These appMWhos, in general, 
deal with the stability, response time, and amplitude determined by solution of deferential equations as 
discussed previously. These solutions fall into the categories classically called homogenous and nonhomo- 
genous and typically are called stability and response, Determination of these characteristics in the 
stability area can be accomplished using (1) Routh criterion, (2) root locus, (3) Nyquist, and (4) Nichols 
techniques. Complex variable theory of functions is the underlying theory used in tiiese techniques, 
such as Laplace transforms. Closed form solutions to the equations can be used in special cases, but are 
not a generally available technique, Numerical integration of the equations is generally applicable and 
with modern computers is efficient for large systems. Analog or hybrid computers also are an efficient 
approach, particularly for nonlinear systems. 


Recognizing that all dynamic models of structural systems and control systems are usually cast 
in differential equation form, the techniques just described arc applicable, The major differences between 
the two disciplines arc in the way these equations arc formulated. 


Control engineers are concerned with several aspects of the structural dynamic interaction prob- 
lem, First, control feedback logic is used as a technique for placing a dynamic system at a given state, 
normally a displacement or rotation. Tills is accomplished using a response command that drives the 
system to tills position by actuating a control force, If the system is stable, then the dynamic system 
arrives at the desired position. Tliis involves a position command about some static or normal position, 
Secondly, in achieving this position command, the engineer is concerned with the rate of reaching this 
position and the overshoot errors and recovery time produced in changing the dynamic system state, 
Thirdly, closely related to this is the stability of the system due to introducing control. Finally, he is 
concerned witli the response of dynamic systems to any external environments that result in perturba- 
tions to the response position he is trying to achieve. These four areas exist whether one is concerned 
with using control to optimize structural design or the control of a system to a desired state which results 
in dynamic responses that arc unwanted and must be contained or reduced (stability augmented systems, 
etc,), When these areas are considered, the basic vibration response characteristics discussed earlier result; 
however, it can be advantageous to formulate the equations in terms of a transformed or different 
coordinate system, Tills coordinate system is formulated in terms of an error coordinate which is the 
difference between the commanded state and the response state, normally called Xj n - X out - e. As a 

result, special analysis tools have been developed around these redefined equations which allows quick 
insights into the control system parameters in terms of the system stability and response. Obviously, one 
does not have to go to this special formulation, but can use the basic structural dynamics formulation 
given previously where control is an input force proportional to state. The latter formulation results, in 
many cases, in a less efficient analysis approach, It conserves, however, all the structural dynamics 
engineer’s built-in understanding and intuition. Therefore, the control engineer needs awareness of this 
data base. At the same time, the structural dynamic engineer must strive to reinterpret his firm data base 
in terms of the transformed equations approach used by control so important to understanding and 
implementing control concepts. Extensive literature exists on these techniques from both disciplines 
and is not a focus of this paper; however, the following example depicts the general concept. 


Using the original system with some different control concepts will allow reformulation of the 
equations. A differentiator will operate on a position signal to obtain a rate in such a manner that the 
resulting equation has the form: 
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e * ~C e e - -C e (X in - X 0 ) 


( 20 ) 




Introducing an electrical amplifier as a means of changing signal, then control force becomes for 
a control amplifier of gain m 

M (a 0 e + aj C c e) F* (21) 

Then, 

m X 0 + C X 0 + K X 0 « f(t) + /xF* (a 0 e + aj C c e) (22) 

which becomes 

MX 0 + (C + C e M F* ai ) X 0 + (K + /iF* a 0 ) X Q = C e (X F* aj X Q + /xF* a 0 ^ + f(t) (23) 

Sotting f(t) = 0 and putting this equation in transform form, the response can be written as. 
(C e ai6 + a 0 )/iF* Xj 

X st 1 z. : (24) 

° M5 2 + (C + C e m F* aj) 5+(K + /xF* a 0 ) 

Or in conventional vibration form.' 


CC eai 5 + a 0 ) MF^/m Xj 
X Q 1 

6^ + 2foo6 + 

where 

C + C e mF* ai 

2 fu = 2 1 

rn 


and 



K + juF* a Q 


(25) 

* 


(26) 


( 27 ) 
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which is the same as equation (20) when the acceleration term is set to 0 and a 0 F* and ajP* *s rede- 
fined to have an amplifier gain, This system is shown in block diagram form in Figure 8. 



Figure 8. Closed loop control augmentation diagram. 


Equation (25) is now in a form that is easily ammenable to various control commands, response 
augmentation, and analysis techniques, achieved through the introduction of a command signal and feed- 
back error signal while maintaining the basic physical system presented earlier. As discussed previously, 
what one is dealing with is the transformation between the time domain and frequency domain or vice 
versa using a transformed parameter e. This formulation, therefore, allows the use of all the techniques 
developed in the field of operational mathematics as well as control theory. It is assumed that the 
reader is knowledgeable of these techniques for (1) generating time responses from equations in operator 
forms using inverse transforms, etc!^(2) generation and understanding of transfer functions Y(S) or 
Y(io>), and (3) numerical integration of differential equations, The area where, in general, the engineers 
lack understanding or insight into control analysis techniques is the stability analysis techniques discussed 
earlier. Whether an engineer conies up through the control or structural dynamics side, he knows that 
the stability of a dynamic system is obtained from the homogenous solution through the roots of the 
characteristics equation. The loss of intuition arises because of the special formulation of the equations 
used in control theory. 

The structural dynamics engineer basically thinks in terms of frequency and damping which is 
preserved only with the root locus technique. Other techniques for stability describe stability margins 
in terms of phase and gain which arc meaningful terms but at the expense of the concrete damping 
expression. As was shown in equations (25), (26), and (27), the resulting transformed equation has the 
form for the transfer function: 


* 

X 


b P m + + b 

0 _ Y (P) = m 0 

a n P 11 + ... + a Q 


(28) 


and the closed-loop differential equation is 

(a n P n + ... a 0 ) b 0 » (b,,, P m + ... b 0 ) 0; 
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Therefore, the homogeneous solution that determines stability is arrived at by setting the right-hand side 
equal to zero and solving for the complex roots 6j’s where 


6 q a ± ico (29) 


The problem is not writing the characteristics equation but finding its roots, The characteristics 
equation can be written in the form of a polynomial, or can be expressed in state form, or as coeffic- 
ients of second order different!? equations in matrix form, The larger the number of degrees of freedom 
or the order of the polynomial, the more difficult the solution for the roots are, This results in the use 
of transformed equations in conjunction with the many techniques based on complex variable theory 
that use the open-loop transfer function as the basic equation formulation, These techniques arrive at the 
answers graphically or numerically in tire frequency domain instead of using polynomial solvers or 
matrix root iteration techniques. Tire power of these techniques such as Routli criterion, Nyquist 
stability criteria, Is that the stability boundaries (positive real roots) can be found without obtaining 
roots; however, as stated previously, tire absolute level of stability Is not known, The Routli criterion Is 
based on the mathematical principle that the coefficients of the characteristic equation are the sums and 
products of the roots of the equation, When carried through to completion, it is easily shown that for 
stability all the coefficients must be positive, For simple characteristics equations (low order, few 
parameters), the coefficients can be written in terms of key parameters and thus the stability boundaries 
can be plotted versus parameter variations. Using equation (24) as an example, the coefficients are? 


a >5 ® m 


a I a C + C c mF* nj 


(30) 


n 0 ® K + juF* a Q 


which says that m must be positive which physically is correct and that: C *1* C e juF* nj and K + juF* a 0 

must be greater than zero, This gives two equations in two unknowns which can be solved as a function 
of the other parameters, Since there are no cross couplings, the solutions are Independent which means 
the damping and frequency and thus the stability can be determined independently, in most eases, the 
system is coupled containing many parameters which must be balanced to achieve the most optimum 
stability characteristics. Literature is full of the application of Routli criterion, 


Nyquist stability approaches are also well documented, It is based on complex variable theory 
which says that if you plot an envelop of points in the complex root plane (closed contour) and calcu- 
late the value of the complex function (transfer function), a closed contour will also be developed in the 
f(6) plane. This closed contour about the origin will circle the origin a total of the number of poles 
minus the number of zeroes where the poles are the roots of the denominator and the zeroes are the 
roots of the numerator of the transfer function, The stability can be determined therefore by graphic 
or numerical techniques and does not require finding the roots of the characteristics equation. In general, 
the equations are transformed such that the transfer function has the form: 
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m 


Y 0 (s) 


x ° * TTy^s) x in 


(31) 


In this case, the number of contours encircled minus one instead of the origin. Root locus tech- 
niques using present day computer technology has not placed as much emphasis on graphic techniques 
as when the engineer had to determine system stability, sensitivity, and margins without high powered 
computer technology and was forced to use graphic techniques. Root locus techniques plot the position 
of poles and zeroes of the closed loop transfer function as a function of some control parameter such as 
the amplifier gain (x , which is a gain in the open loop transfer function. In other words, the root locus 
procedure starts witli the open-loop transfer function and ends with the poles and zeroes of the closed- 
loop system. Before proceeding it is important to remember that the closed-loop transfer function using 
inverse transforms describes the transient response characteristics. In general form: 


X 0 (t) = Cj c s > 1 + C 2 e Slt + ... C n e s>t 


(32) 


where the poles S u <= o n + i co n determines the stability and frequency of response and the zeroes fix 

the size of the transient term for a particular input, i.e., the constants C|, C 2 , and C n , In this case, 
equation (34) is rewritten: 


X r 


X, 


IC 


G 


X{ 1 + X Q 1 + Kq 


(33) 


where K is an arbitrary gain of the control amplifier. 


The characteristic equation becomes 


(P - Sj) (P - s 2 ) ... (P - s n ) = 1 + X 0 (S) = 1 + G(S)K = 0 


(34) 


Graphically, then one plots the points where 


G(S) = ~ 


(35) 


Using complex variable theory, the solution to equation (35) is found where the angle of G(S) 


Angle of G(S) = angle of - ■— = it + K 2tt 

K 


( 36 ) 
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(37) 


K - integer and the magnitude of G(P) 


K 


Thus, the angle of G(S) 

2 P 

Angle G(S) = <A 4 £ zero vectors - £ pole vectors 


(38) 


Magnitude G(S) » |A| 


Product of Z zero vector magnitudes 
Product of P pole vector magnitudes ’ 


where these equations are derived from the factoring of 


G(S) - 


A(S-Z 1 )... (S-Z n ) 

cs - Pji' .„'cs - p n )' 


(39) 


Therefore, using complex variable theory, one can treat these complex members as vectors and 
graphically evaluate G(S), It is not the purpose of this memo to go further into the details of root 
locus, since it is well documented. Modem computers by-pass this graphics approach and solve for the 
roots numerically which is a more efficient way. It is clear the solutions obtained in this manner are 
equal to and consistent with those normally dealt with in structural dynamics. As stated previously, 
modem day approachs do not depend on these graphic techniques; however, it is imperative that the 
underlying principles inherent in these approaches be understood as well as staying abreast of the data 
base developed in conventional structural dynamics analysis. 

Control theory is not ideal in application. Forces cannot be generated without introducing phase 
lags. Also, through the use of electrical networks, and a digital control theory, the control engineer can 
generate a control force or control signal that changes in amplitude and phase as a function of frequency 
given additional flexibilities and power to design engineers. These effects do not change the basic con- 
cepts presented. They do increase significantly the complexity of obtaining the solutions. 

Section IV looks at these concepts in terms of the response of an elastic space vehicle or aircraft 
response to environmental excitation. 


SECTION IV. SPACE VEHICLE ELASTIC BODY RESPONSE CONCEPTS 


Before going to elastic body response, it should be mentioned that the rigid body response of a 
space vehicle in rotation and translation are mere extensions of the concepts just presented where the 
major stiffness and damping are due to the control system with some augmentation from aerodynamics. 
The major mass moment of inertia portion comes from the structural configuration with little or no aug- 
mentation from control. This subject, rigid body response, has been treated extensively for many years 
[ 1 ] and is not repeated here. This section will discuss only the elastic body effects because of their 
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strong interaction with control and performance issues of large space structures, where not only control 
of response but control of curvature or shape is important. In these cases, requirements on ferrying 
weight and overall size result in a very light, large, and low frequency structural system. These 
constraints place strong emphasis upon a coupled system optimized design approach. Also, during ascent 
flight and aircraft flight, large aerodynamic and gust loads and uncomfortable ride qualities result from 
the elastic body transient response to these disturbances, Use of control to reduce these responses is a 
key design area. 

As all structural dynamicists are aware, the solution of the structural dynamic equation with no 
external forces leads to a set of normal modes and frequencies commonly called eigenvectors and eigen- 
values, The modes are orthogonal and thus uncoupled. Since each mode can be characterized as an 
effective mass and stiffness, called generalized, one can then couple them with external forces, including 
aerodynamic and control system using energy approaches such as LaGrange’s equations. Using these 
approaches, including generalized coordinates, the equations are in the same general form used in vibra- 
tion theory of both the single degrec-of-frecdom forced response and multidegree-of-freedom forced 
response systems. Extensive documentation is available depicting these approaches and solutions used. 

To understand clastic body responses (loads and stability), the assumption will bo made (later 
removed) that one elastic body mode is uncoupled from the other and that the rigid body angle of 
attack and engine deflection act as known (time wise) forcing functions to this model (Fig, 9), 



In order to write the equations under these assumptions, the gitnbal engine generated force is split 
into two parts; the rigid body generated control induced force is a known function of time, and the 
elastic mode introduced control force is a function of the elastic body modal state (deflection, rate, and 
acceleration). Phasing between the rigid body generated engine and aerodynamic forces is neglected for 
simplicity. When phasing is neglected, the equation for a bending mode is written as follows: 
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where 


r? M Ct) = bending mode generalized coordinate, 
co Bm = bending mode natural frequency, 

M b = bendinf mode generalized mass, 
f B;U - structural damping, 

" local angle of attack aerodynamic term, 
= local angle of attack aerodynamic term, 
53 rigid body aerodynamic force term, 

F s = veliicle tlirust, 

Yj 7 = mode deflection at engine. 


( 40 ) 


(41) 


To the uninitiated, the engine nozzle is gimbaled to produce a lateral control force proportional 
to the tlirust times the sine of the deflected angles. For small angles, this force is equal to the tlirust 
times the angle in radians. Assuming that ^elastic rcsu Es from signals arising from body-fixed accelerome- 
ters, rate gyros, and position gyros, the control equation becomes: 

Elastic * a 0V (t) Y ' (X g ) + a l V° Y 'C X R) + *2* COY (X A ) , (42) 
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where hq is position signal gain, a\ rate signal gain, &2 accelerometer signal gain, Y'(X g ) the bending 

mode slope at the position sensor, Y'(Xr) the bending mode slope at the rate sensor, and Y(X^) the 
bending mode deflection at the acceleration sensors. 

Substituting equations (41) and (42) and simplifying result* gives: 



m B - S - “1 Fs y E y '(Xr>1 V‘> 

+ Mb - »2 P s Y E y < x A> 


[<V M b - B„ - a 0 F s Y E Y'(Xg)] r,„(t) 
+ «B - “2 F S Y E y < x a> 


Q(t) 

[Mb - «2 f s y e y < x a» 


(43) 


It is clear that the above generalizations were made for one sensor of each type; however, the use 
of more than one sensor does not destroy the use of the analogy, since the total signal is the sum of the 
voltage coming from each control loop. The effects of multisensors on the roots, and therefore the 
cross-coupling between modes, etc., will be addressed later. How the response is altered by the control 
system on the roots is now discussed for the ideal case. The equation becomes: 


(t) + 2S li(t) + (« 2 + (3 2 ) (t) = R Q(t) 


(44) 


where 
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kr (Apparent Damping) 


(46) 


(or + j3 2 ) = [co B 2 M b - B m - a 0 F s Y E Y'(X g )l R 2 = w Q 2 (Apparent Frequency) . (47) 

M 


The roots to the equations are obviously <x ± p defined in equations (46) and (47). Using these 
roots, the solution to equation (44) gives insight into response-relieving mechanisms. Typical solutions 
can be found for constant coefficients and known Q(t)’s. Assuming that Q(t) is some known (Laplace) 
transform, then 
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for initial conditions equal to 0. Letting Q(t) be a ramp input or Q(t) “ K 3 t, then 
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Finally, if Q(t) is a sine function 
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7,2 
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e“^ sin fft + ^p)3 


where 


1//9 ~ ^7 - tan" 
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a 2 -? 2 


(56) 


( 


As expected, all cases show a difference between acceleration and displacement of the co 2 factor 
on the transient part of the solution except for the sinusoidal forcing function which also contains a 

steady-state term with a factor of S2* difference. Considering the solution to the ramp, step, or impulse, 

the magnitude of the constant can be changed by the term R 2 by use of accelerometers. These solu- 
tions — ramp, step, and impulse - can also be altered through each of the sensors as they alter the 
frequency or damping [equations (46) and (47)3. Rate gyros change the damping of the system either 
positively or negatively depending on the sign of the modal deflection Yg and Y'(Xj^) and the rate gyro 


gain, aj. Choosing the sensor location or gain such that a increases results in greater damping and lower 

transient. Choosing an accelerometer location and gain such that R increases, increases both the damping 
and the frequency, thus allowing the accelerometer to be used as a modal suppressor from both the 
damping term and the frequency term, Position gyros can be used to alter the frequency by a proper 
choice of the sensor location or feedback gain, ap. The amplitude response (both steady-state and 


transient) in these cases is reduced if the frequency is increased; however, the accelerometer output is 
proportional to the frequency squared times the transient portion of the solution. Increasing the 
damping lowers the peak transient response. All three types of input forces are expected during flight 
since the wind contains some form of each type of input. The response (acceleration or amplitude) can 
be reduced by increasing the frequency or damping. 


A more important type of force from the bending mode standpoint is the sinusoidal input. This 
represents the turbulence portion of the atmosphere, which can have frequency content in resonance 
with the bending mode. Also, not only is the transient term important but the steady-state term can be 
of a larger magnitude in both acceleration and amplitude. Again, increasing the damping decreases the 
amplitude and thus reduces the transient response. Increasing the frequency may not be feasible; how- 
ever, since the resonance term contains^ 2 + a 2 and £2 2 , and the amplitude of the frequency increases as 

j3 c + or and approach equal values. In this case,^the frequency shift must be chosen to detune the 
system from the forcing frequency. Additionally, for this case, the accelerometer can be used to reduce 

the overall amplitude through R 2 which multiplies the solution (Fig. 10). Care must be exercised in 
using tills term for reducing amplitude when at the same time it may increase the amplitude through 
either finer tuning (with forcing function), or decreased damping and frequency* The change in damping 
and frequency can be obtained by using the various sensors as discussed previously. 

The discussion thus far has indicated that accelerometers, rate gyros, or position sensors, can be 
used for augmenting elastic system response. The results point out very clearly that one must have a 
very accurate description of both the vehicle modal characteristics and the input force (wind) to effec- 
tively design a control system for modal suppression using these types of sensors. Also, the coupling of 
these sensors in the rigid body control (used as input force for mode) is very important and cannot be 
neglected. 
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Figure 10, Modal suppression. 

The previous interpretations can be stated in another way. The basic notion here is the freedom 
offered by a sensor complement in locating closed-loop eigenvalues as a possible source of quality 
measures. Tills is motivated by two considerations. First, classical experience with root loci and fre- 
quency domain design techniques provides tested insightful relationships between the performance cap- 
abilities of a controlled system and the closed-loop pole arrangements permitted by sensors, Such notions 
as stability, frequencies of oscillation, damping of individual modes of response, and dominance are all 
apparent from the pole constellation. Secondly, there is a fundamental connection between pole place- 
ment and the concept of controllability. 

The previous discussion was based on the assumption that ideal control signals and response exist. 
Although tins is not true, the principles remain the same as long as the gain and phase lag changes that 
take place in reality are considered. Also, the assumption is made that each mode is completely inde- 
pendent of the other, wliicli is not true. To illustrate tills, a two-sensor case will be presented first, then 
a two-mode case, 

If two accelerometers are used instead of one, the denominator in equation (47) becomes: 


= M B - a 21 F S Y E Y(X A] ) - a 22 F s Y E YQC^) 
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This allows a choice of gains anti sensor locations that would cancel the accelerometer effect or allow 
any mixture of effects (gains) between accelerometer locations, The other coefficients in equation (47) 
could be modified in the same manner by using two or more rate or position gyros, This not only 
illustrates the complexity of using many sensors but also the flexibility. 

Extending the concept to two bending modes, but neglecting certain rigid body coupling, results 
in the following equations which are derived by assuming only one sensor of each type in the control 
equation, The control equation (elastic body feedback portion) is 


^elastic = *0 «i y i'< x r) + fcY, '«,*>] 


I 


+ g 2 « 1 Yi<x a ) + %y 2 cx a )] . 

The coupled bending dynamics equations given in matrix form, using this control law are as follows: 


(58) 


* 


" m B1 -*2 F s Y E1 y 1 < x A> 

- 62 F s Y E1 y 2< x A> ’ ( 

04 
* 1 

- 82 F s y E2 M < x A> 

m B 2 - 8 2 F s Y E2 y 2 < X A> j 

h 

, 

2r B i M b) - C, , - aj F, Ygj Y | f (X R ) 

*" a i F s Y E2 Y 2 f ( X R) " c 12 



“ a l F s Y E2 Y l'C x R> - C 21 

2 ^B2 w B2 m B 2 “ c 22 " a l F s Y E2' X R_ 


^2 

M B) - Bl 2 - B„ - o 0 F s Y E , Y,’(X 6 ) 

- B 12 - a 0 F s Y El Y 2 < X g) 


*1 

- U 12 - a o F s Y E2 Y l'( x s> 

m B 2 "B2 2 - B 22 - a 0 F s Y E2 Y 2' < X g> _ 


*2 



(59) 


The coefficients show that what is done with one sensor for one mode can be offset by the 
redundant signal from the second mode. If the system is extended to include many sensor gains and 
force input locations, the tradeoffs arc apparent but too difficult to illustrate. Although the concepts 
for one mode hold for this more general case of two modes, the design problem is increased many fold 
because of cross-coupling through the control system. Obviously, the tilings that help suppress one mode 
could easily aggravate another, With several modes and sensors, a procedure must be used that provides 
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insight into Important characteristics and that gives first cuts at the gains and sensor values and locations, 
Numerical response techniques, optimal control techniques, sensitivity studies, and root locus techniques 
(eigenvectors and eigenvalue routines) allow the development of trades between control force locations, 
sensor locations, and control law logic that provide the desired response in each mode or the over 
response as a deflection and load, 

Illustrated in this section has been the extension of the concepts discussed previously to elastic 
body response of a space vehicle. Additional concepts of modal gains and multisensor control logic as 
well as extension to multimodcs systems were addressed, all prescribing the response in the general form 
a structural dynamleist is familiar with, Non-ideal control was not addressed; however, It is easily 
included by writing the differential equations for these effects and adding them into the matrices and 
generalized coordinates given in equation (30) [2], 


SECTION V. NON-IDEAL CONTROL EFFECTS ON THE SOLUTION 


As stated previously, the control system docs not perform in tliis manner, Also the control 
engineer can alter the gain and phase electronically, producing greater flexibility to alter structural 
response, Therefore, the analog presented must be extended to cover these conditions, Assumptions 
must be made in this case also; namely, that these changes can be made without introducing excessive 
degrees of freedom and that those introduced do not create unstable modes at some other frequency, 
all of which can be accomplished with proper design attention. Also, It is assumed that the coupled 
frequency is fairly close to the uncoupled frequency of the spring mass system, This assumption can be 
removed later, Using these assumptions and going to the equations using the control equation in the 
form 


( 60 ) 


where Kj is the gain and 0| the phase of the control system near the natural frequency of mass spri g 

system cither due to inherent lags or is artificially introduced to achieve a desired response. For each of 
these terms, they can be rewritten in complex form m 


Kj (cos 0 (oc) + i sin 0j (cc)) 


s 

(61) 


It is desirable to remove the complex term (i). This can be accomplished under the assumption 
already states that the frequency ./ill be near the natural frequency by dividing by co and taking a time 
derivative of the variable. Performing these operations and assuming that the acceleration can be approxi- 
mated as -co^X, the control equation becomes 


7 - F* ag ICj [X cos 0j(co) + X/co sin 0j(w)] + F H ' aj [X cos 02(co) + X/co sin 02(co)] 


+ F* & cos 03 “ X/co 03 ( w )l 
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Reinserting these terms back into equation (42) gives 
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Tire power of control becomes obvious when it is recognized that the K’s and 0’s can be determined 
using state-of-art filter-design techniques, The <^’s can take on any value from 0 to 2 tt, the K’s are usually 

values between 0 and 1. The choice of the signs and magnitudes on the control gains, aQ, aj, a 2, coupled 

with the choices of K’s and 0’s allows altering the response of the basic equations in many ways: (1) by 
changing the structural dynamic design, (2) through basic control system gain choice, and (3) by intro- 
ducing filters or networks that change the feedback as a function of frequency. In all cases, the basic 
understanding reached for the single degree of freedom system is preserved. Obviously, one cannot get 
these frequency dependent changes without altering the number of degrees of freedom; however, if 
proper attention is given, it will not destroy the analog given. 


In summary, up to tills point in the report, one sees that a control system in all its basic linear 
forms can be interpreted as pseudo structural elements, preserving the general response characteristics 
which a structural dynamics engineer is so familiar with. Ideally, then control can be interpreted as a 
means of adding or subtracting stiffness, damping, and facatis, all gained without changing the structural 
design. Adding the influence of filtering networks, any feedback term can be made to act as springs, 
dampers, or mass. This opens up a whole new vista of design approaches that includes substituting 
control for structure and is the basis for the relative new field of control configured vehicles. The 
penalty paid for these additional design parameters is the redundancy required in the control system for 
fail operational or fail safe requirements. 


SECTION VI. ANALOGIES TO VARIOUS CONTROL CONCEPTS 


So far, all the discussion has centered on the use of linear, proportional, central, control theory. 
The analogies shown do not preclude, however, nonlinear control or decentralized control, Here, the 
structural dynamicist would just move into nonlinear vibration theory or specialized analysis and develop 
the same kind of analogies. 

Digital control theory can usually be put in the form of the continuous analog systems used so 
far in this paper if the sample rate is high enough relative to the structural frequencies of concern. This 
is always a good starting place. From that point on, new analogies need development but fall into the 
categories of sample data system analysis techniques currently employed and sample data control theory 
analysis techniques, which through transformations put the operations back in this same general form 
discussed previously. Certainly, one can think easily in terms of RMS, spectrums, etc., and how these 
are altered for dynamic system with control. The development of various analogies are left to the reader. 

Force point sensing is an old concept; however, it has new emphasis in that by using distributed 
control (at the force inputs) the concept can be implemented for numy modes. The idea is to sense and 
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control the response at the same place on the structure. Ideally, tills system is always response reducing 
or stabilizing; however, control force actuators introduce phase lags that can change this effect, Intro- 
ducing a filter or network that compensates for the actuator characteristics would put one hack to the 
ideal case discussed. 


The concept of distributed control is to put the logic, sensors, and force actuators in selected 
areas of the structure. This means the logic and desired response for each area must be well understood; 
otherwise, cross talk between each area creates real problems as they compete and couple with each 
other. This approach does not, however, destroy the analogies developed so far, This concept Is very 
powerful for shape control. Using the bending mode equation (14) derived in the previous section 
requires only that terms are included for each sensor, control force, and logic accounting for the slopes 
and deflection where each element is located. The only complexity is the additional terms and ability 
to differentiate their Influence, The problems get very messy when going to a multimode system; how- 
ever, many more variables are available to got the optimum response solution, 


Putting the centralized control and distributed control together and allowing the centralized 
control to be hierarchal gives additional flexibility but at the expense of complexity (Fig. 11). The 
reader can easily develop these analogies if desired. 


□ CONTROL SENSOR 
O CONTROL EFFECTOR 



Figure 11, Distributed control system. 




SECTION VII, PROBLEMS, CONCERNS, AND FUTURE CHALLENGES 


Everything discussed has made the tacit assumption that the sensed state is ideally identified, 

This is not the case in the real world. Many times noise, other subsystem response frequencies, etc., 
obscure the desired signal, The whole area of state identification and pattern recognition must be 
extended. With bolter and better sensors, more complex control logic, and lvigli structural modal density 
of very complex and highly coupled modes makes this a prime area to attack. 

Optimal control approaches have basically dealt with the system response with no weighting 
factors for structural design changes. These areas need to be integrated further than has been accomp- 
lished to date and development of a means of doing total system optimized design in terms of structure, 
materials, thermal, and control. 

A major problem facing tliis new discipline is the basic sensitivities of the responses, design, etc,, 
to the uncertainties of both the structural model and the control system model, Here is needed learning 
control devices, desensitizing approaches, and sensitivity analysis techniques for large parameter systems. 

Techniques for truncating or selecting critical structural modes and control system characteristics 
is one of tho major problems facing this new discipline, If the problem solutions are reduced to a 
manageable size and complexity, major effort needs to be started in both criteria and reduction tech- 
niques developed. 

Obviously, in the future not only space vehicles and large space systems but the whole field of 
energy and transportation must deal with this new discipline. The future challenge is to either create a 
now discipline or cause the separate disciplines to understand, communicate, and work together in a way 
not generally done in the past. It is hoped that tliis short introduction to some of the analogies will 
cause someone to look afresh at these areas. 
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